function A = plotrho(stop)
for i=1:1:(stop-1)
    [filename,err] = sprintf('rho_%d.dat', i);
    A = load(filename);
    plot(A);
    %ylim([-1e-6 0]);
    pause(0.01);
    %hold on
end
